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We analyze the topological transition and localization evolution of disordered two dimensional 
systems with non trivial topology based on bipartite lattices. Chern insulators with broken time 
reversal symmetry show non standard behavior for disorder realizations selectively distributed on 
only one of the sublattices. The Chern number survives to a much stronger disorder strength (one 
order of magnitude higher) than in the equally distributed disordered case and the final state in the 
strongly disordered case is metallic. 


I. INTRODUCTION-SUMMARY 

The metal insulator transition and the associated is¬ 
sue of the existence or not of metals in two dimensions 1 
have fascinated physicists in the last century and have 
been the focus of intense debates and great advances in 
the understanding of low dimensional material physics. 
The interest of these materials is of special relevance 
today due to the unusual applications expected follow¬ 
ing the experimental capability of synthesizing and ma¬ 
nipulating layered compounds. 2 But also the fundamen¬ 
tal physical aspects are still providing surprises mostly 
around the Dirac physics based materials and their re¬ 
lated topologically non-trivial features. The issue of how 
electrons in a metal become localized by disorder in the 
presence of non-trivial topology started with the anal¬ 
ysis of Landau levels in the integer quantum Hall effect 
and was rapidly adapted to other time reversal preserv¬ 
ing topological matter. 3 ® The standard classification of 
disordered classed® based on the Alt land-Zirnbauer sets 
of random matricel 9 was completed to include topolog¬ 
ical features and the “tenfold way” was set in Refs. ITOl 
and [IT, 

Two dimensional (2D) topological insulators belong 
to class A (unitary) with time reversal symmetry bro¬ 
ken or All (symplectic) time reversal invariant. Promi¬ 
nent examples in these classes are the Haldane 12 and the 
Kane-Mele model! 13 They are characterized by a quan¬ 
tized Chern (or spin Chern) number C. The localization 
transition in these cases is accompanied by a topological 
transition between a topological (band) insulator and a 
trivial (Anderson) insulator. The topological transition 
of these classes is assumed to be well understood. The 
standard mechanism is referred to as “levitation and anni¬ 
hilation” of extanded statesP® For moderate disorder, an 
impurity band of localized state forms inside the gap and 
simultaneously the high energy states in the edges of the 
conduction and valence bands start to localize. As dis¬ 
order increases, the gap is totally populated by localized 
states and the bulk extended states above and below the 
Fermi level carrying the Chern number shift toward one 
another and annihilate leading to the topological phase 
transition. The transition from a Chern to an Anderson 


insulator under a random potential disorder has been an¬ 
alyzed in Ref. [15. In the standard classification there 
are also the chiral classes characterized by having sub¬ 
lattice symmetry. They are typically a naly zed by adding 
random disorder to a bipartite latticed 6 17 Graphene be¬ 
longs to this class and its localization properties have 
been studied at length.^^I The localization behavior of 
Chiral classes is more complicated but they do not have 
associated topological invariants in two dimensions. 

In this work we focus on the sublattice structure of 
the honeycomb lattice underlying 2D topological insula¬ 
tors and show that the critical disorder which signals the 
appearance of the final regime has an extraordinary sensi¬ 
tivity to the disorder selectivity over sublattices/orbitals. 
The Anderson transition in the chiral classes is particu- 
larly interesting with interactions playing a role 17 but, as 
mentioned, there are no topological compounds in this 
class in 2D. We see that, for disorder realizations affect¬ 
ing only one of the sublattices, the Chern number sur¬ 
vives to a much stronger disorder strength (one order of 
magnitude higher) than in the equally disordered case 
and the localization transition never takes place. We ex¬ 
emplify our findings with the Haldane modeP 2 with two 
different types of disorder: Anderson disorder and va¬ 
cancies. It is interesting to note that, although there are 
qualitative differences on the spectral evolution as dis¬ 
order is increased, the two cases behave similarly and 
the important difference lies in the two possible disorder 
realizations. 

We visualize the topological transition as a transition 
from the honeycomb (topologically non-trivial for finite 
Haldane mass) to the triangular (trivial) lattice within 
the same symmetry class. The “localization transition” 
is more complicated but in the limiting cases it occurs 
between a band insulator (topologically non trivial) and 
a metal, again without changing symmetry or dimension¬ 
ality. 

II. MODEL AND METHODS 

We use the Haldane model 12 ! as a generic example of 
a topologically non trivial system in symmetry class A. 
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It was initially proposed as a tight binding model in the 
honeycomb lattice with complex next nearest neighbor 
(NNN) hoppings carefully adjusted so as to mimic a mag¬ 
netic flux flowing in opposite directions through the two 
triangular sublattices. It was build to show the possi¬ 
bility to have Hall conductivity in zero external mag¬ 
netic field and it is nowadays the prototype of a modern 
topological insulator. Interestingly, a physical realization 
of the model has been implemented in Ref. [22] (see also 
Ref.;23). The Haldane model tight binding Hamiltonian 
can be written as 

h = ~tj2 4°j + -*2 e ~ l4>ii 4 c i (!) 

{id) {{id)) 

+ M ^2 rjicjcj + H.c., 

i 

where q = a, b is defined in the two triangular sublattices 
A, B that form the honeycomb lattice. The first term t 
represents a standard real nearest neighbor hopping that 
links the two triangular sublattices. The next term rep¬ 
resents a complex NNN hopping acting within 

each triangular sublattice with a phase 4>ij that has op¬ 
posite signs (j>ij = ±0. This term breaks time-reversal 
symmetry and opens a non-trivial topological gap at the 
Dirac points. We have done our calculations for the sim¬ 
plest case 0 = 7r/2. The last term represents a staggered 
potential^ = ±1). It breaks inversion symmetry and 
opens a trivial gap at the Dirac points. 

We will discuss two types of disorder: Anderson and 
vacancy disorder equally or selectively distributed among 
the two sublattices. Potential (Anderson) disorder is 
implemented by adding to the Hamiltonian the term 
^2 ieA b e i c l c ii with a uniform distribution of random lo¬ 
cal energies, Si G [—W/2,W/2\. For selective disorder 
the sum runs only over one sublattice. The same is done 
for vacancies which are introduced by removing a given 
lattice site and the hopping to its neighbors. 

The Haldane model belongs to symmetry class A where 
the different topological phases can be characterized by 
a Z-topological number, the Chern number c! 10 In the 
clean system it can be computed from the single particle 
Bloch states u n ( k) as: 

Cn = Y f ( 2 ) 

J S 

where the integral is over the unit cell S and fTJ(k) 
is the z component of the Berry curvature, U n (k) = 
V k A*4 n (k) defined from the Berry connection: *4 n (k) = 
(iz n (k)| — iVk|^n(k)). When translational invariance is 
broken by disorder or any other perturbation, the Chern 
number has to be computed numerically in real space. 
The problem has been addressed in the literature and a 
number very efficient numerical techniques based on the 
use of twisted boundary conditions are now available!— 1 — 1 
From the computational point of view an efficient ap¬ 
proach was developed by Fukui et al. 25 where Chern num¬ 
bers are calculated on the basis of a discretized Brillouin 


zone. Based in this technique Zhang el al. 24 have recently 
proposed an efficient method applying a coupling-matrix 
approach to calculate the Chern number in finite disor¬ 
dered systems. The method is faster than previous ones 
allowing to obtain the Chern number with a single diag- 
onalization for the Ni x AT 2 system, if TVi, Ah are large 
enough. The technique has been probed to give very 
good results for the Haldane model. In the present work 
we have calculated the Chern number for the bulk spec¬ 
trum following the coupling-matrix approach described 
in Ref. [24l 

The bulk density of states (DOS) have been computed 
using the recursion Green’s function method as used for 
example in Refs. [26] and [23 


III. RESULTS 

As shown by many authors, a single vacancy induces 
a midgap bound state in the Haldane or Kane Mele 
model,® 33 while a density of vacancies induces an im¬ 
purity band inside the gap. 33 Increasing the density of 
vacancies enhances the spectral weight inside the gap, 
which eventually closes for high enough vacancy density. 
In contrast, potential disorder localizes states near the 
band edges and generate a mobility edge.^ Lifshitz tails 
eventually overlap with incresing disorder and the gap 
closes. We see that Anderson disorder (local random po¬ 
tential) and vacancies show quantitative differences only. 
The crucial difference is found between the two disorder 
realizations analyzed: selective and equally distributed in 
the two sublattices. We summarize our results in Figs.[l]- 
[ 3 ] The evolution of the Chern number with disorder 
strength encodes the topological transition. The vari¬ 
ous cases are shown in Fig. [l] The bulk DOS contains 
information on the spectral gap and is displayed in Figs. 
[2] and [3] 

We will first discuss the situation when impurities or 
vacancies are equally distributed in the two sublattices. 
In Fig. [lja) we show the behavior of the Chern number 
with potential disorder as a function of the linear lat¬ 
tice size d for increasing disorder strengths W in units 
of t. We see that the Chern number is already well de¬ 
fined and stable at very small lattice sizes d. As shown 
in Fig. [jjb), where we plot the Chern numbei^ for the 
largest simulated lattice size d as a function of disorder 
strength W, the obained number ceases to be quantized 
at a mild disorder strength (At < W c < 5£), a behav¬ 
ior already found in the literature ! 24 * 36 * A comparison of 
the DOS of the Haldane model with Anderson disorder is 
shown in Fig. [2ja) for equal disorder and in Fig.[2|b) for 
selective disorder. The three curves shown in each panel 
refer to three disorder strengths W/t = 2,4,6. For weak 
disorder strength W = 2t (blue), well below the critical 
value for the topological transition of the case of equally 
distributed disorder, the two images are very similar. We 
see a well defined gap at zero energy and two peaks of 
similar intensity in the conduction and valence bands, 
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Figure 1. Evolution of the Chern number for Anderson 
disorder [(a), (b)] and vacancy disorder [(c), (d)]. (a) Chern 
number as a function of the linear lattice size d for disorder re¬ 
alizations equally distributed among the two sublattices. The 
different curves correspond to different disorder strengths W 
in units of t. The error bars stand for the standard deviation 
due to disorder average, (b) A comparison of the behavior of 
the Chern number as a function of the disorder strength for 
the two types of disorder: equally distributed (red squares) 
and selective disorder (blue circles). The inset shows the - 
much higher - values of disorder strength for which the Chern 
number ceases to be quantized, (c) Chern number for equally 
diluted sublattices, (d) Selected dilution. 


a situation compatible with the standard picture of the 
mobility edge. For a disorder strength W = 4£ (red), 
close to the critical value for the topological transition of 
the case of equally distributed disorder, the DOS shows 
a clear difference between the two disorder distributions. 
For the selectively distributed case in Fig.[2](b) the struc¬ 
ture of the DOS is similar to the weak disorder case while 
in the equally distributed case in Fig. |2|a) the gap has 
disappeared and the DOS has flattened. Eventually the 
gap also closes for selective disorder by increasing the 
disorder strength, as already seen for W = 6t (green) in 
Fig. [2](b) . This “late” gap closing correlates with a much 
robust topological phase, as we explain below. 

The behavior of vacancy disorder equally distributed 
among the two sublattices is similar to that of the An¬ 
derson disorder. The evolution of the Chern number as 
a function of linear lattice size d for various densities of 
vacancies nv is shown in Fig. [ljc) . Notice the units in 
the vertical axis. The quantization disappears for a crit¬ 
ical density nv c of around 4%-6% compatible with the 
standard mechanism of “levitation and pair annihilation” 
of the extended states carrying the Chern number. 14 A 
similar behavior has been reported in Ref. [37] for the case 
of vacancies. This mechanism is also supported by the 
analysis of the DOS depicted in Fig. 3](a) for various den¬ 
sity of vacancies randomly but equally distributed in the 



Figure 2. DOS of the Haldane model with Anderson disorder 
with equal (a) and selective (b) disorder distribution. The 
three curves refer to different values of the disorder strength: 
W/t — 2,4,6. The arrows indicate the effect of increasing 
disorder. 


two sublattices. It is apparent that vacancies induce a 
finite spectral weight inside the gap. As the disorder den¬ 
sity increases the impurity band widens and the center of 
the conduction and valence band that host the extended 
states carrying the Chern number (1 and -1 respectively) 
move closer. At a density of vacancies around nv ~ 3% 
the gap closes (inset of Fig. |3ja)). The Chern number 
remains quantized up to a density nv c (4% < nv c < 6%) 
where the “annihilation” occurs. Similar behavior has 
been found in Ref. 37| although the critical density re¬ 
ported there is a bit different. 

We will now discuss the results obtained for disorder 
selectively located in one of the sublattices. A compar¬ 
ison of the behavior of the Chern number as a function 
of Anderson disorder strength for the two types of disor¬ 
der, selective and equally distributed in the sublattices, 
is shown in Fig. [ljb). The transition described above in 
the evenly distributed case (red squares) occurs sharply 
around W c « 5 1. For this value, the Chern number in 
the selective disorder case (blue circles) is still quantized 
to one. The topological transition occurs for disorder 
strengths one order of magnitude bigger (inset of Figure) 
(40 1 < W c < 50£) than in the equally distributed case. 

Similarly, for the case of disorder in the form of va¬ 
cancies selectively distributed shown in Fig. 0d), the 
Chern number is still converging to one for a density of 
vacancies of around 30%. This is to be compared with 
the equally distributed disorder in the two sublattices 
(Fig.0 c)) where the graphic departs from 1 already at 
a critical density around 4%-6%. The DOS for selective 
vacancy disorder, shown in Fig.J^b), is also very different 
than in the previous case. Valence and conduction peaks 
do not move substantially as disorder increases and the 
behavior of the gap between the valence (conduction) and 
midgap bands is also different. As can be seen in the in¬ 
set of Fig. §b) the gap decreases and reopens well before 
the Chern number ceases to be quantized what happens 
at a critical density of vacancies of nv c > 30% (Fig. [ljd)) . 
At this density there is a noticeable gap separating the 
midgap band from the conduction and valence bands. 

The case of disorder distributed in only one sublattice 
requires an alternative explanation. We believe that both 
the topological transition and the localization mechanism 
are different in this case. The robustness to disorder of 
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Figure 3. DOS of the Haldane model for vacancies equally 
distributed among the two sublattices (a) and on one sub¬ 
lattice only (b). The insets show the behavior of the gap as 
a function of the density of vacancies for the two situations. 
For equally distributed vacancies the gap closes for a critical 
density nv ~ 3%, as can be seen in the zoomed in region. 
For selected dilution, the gap closes at the critical density 
nv ~ 12% but reopens again for nv > 12%. Arrows indicate 
the effect of increasing disorder. 


the extended states carrying the Chern number is related 
to the fact that one of the sublattices remains perfect. 
The component of the wave function corresponding to 
the “clean” sublattice stays delocalized in a way similar 
to the vacancy states of bilayer graphene described in 
Ref. [26| The fact that the topological transition in the 
vacancies case occurs for nv c > 30%, while the perco¬ 
lation transition of the triangular lattice takes place at 
50 %JH39 points to the possibility that the topological 
transition occurs when the disordered sublattice is de¬ 
stroyed and the extended (Bloch) state is driven by the 
underlying triangular lattice that is topologically trivial. 

The case of vacancies is very intuitive: In the limit 
when we have removed all the atoms in one of the tri¬ 
angular sublattices what remains is the extended wave 
function corresponding to the triangular sublattice. Let 


us emphasize that we should not confuse the triangu¬ 
lar lattice with homogeneous complex hoppings, which 
break time reversal symmetry, with the familiar triangu¬ 
lar (or hexagonal) lattice with real hoppings. A discus¬ 
sion of the transition between the honeycomb and two 
independent triangular lattices is done in Appendix A] 
The corresponding bands are shown in Figs. [4] and [5 of 
the Appendix [Aj 


IV. DISCUSSION AND OPEN ISSUES 

It is by now clear that topological features can change 
the localization behavior of the disordered electronic lat¬ 
tice systems. The topological transition (topological in¬ 
dex going to zero or to a non-quantized value) and the 
localization transitions are related but do not occur, in 
general, at the same critical value of the disorder. In the 
“ten-fold” way classification of topological insulators,® 
chiral classes (with well defined sublattice symmetry) do 
not support topological features in 2D and the issue of 
selective disorder (affecting only one sublattice) has only 
started to be explored. 21 The topological transition in 
the 2D (non chiral) classes is assumed to occur through 
“levitation and annihilation” of the extended states car¬ 
rying the topological index.® The mechanism is less un¬ 
derstood in the chiral (sublattice symmetry) classes.® 
In the one dimensional chiral case analyzed recently in 
Ref. ST] the topological index C remains quantized and 
non-fluctuating even when the bulk energy-spectrum is 
completely localized and after the insulating gap has 
closed. C changes abruptly to a trivial value when dis¬ 
order is further increased. 

We observe a similar behavior for the Chern number 
C for selective vacancy disorder: C remains quantized 
up to a large disorder strength but in our case extended 
states remain present possibly all the way until one of 
the sublattices is totally destroyed. This type of disor¬ 
der realizes the transition from honeycomb lattice physics 
with topologically nontrivial features to the two topolog¬ 
ically trivial uncoupled triangular lattices. Our results 
for the selective disorder do not convey with the stan¬ 
dard description of the topological transition based on 
“levitation and annihilation” of the extended states car¬ 
rying the topological index. 

Sublattice symmetry plays a very subtle role in 2D 
topological insulators. It is a necessary ingredient in a 
minimal model to support non trivial topology (a single 
band is necessarily trivial) but at the same time non triv¬ 
ial topology requires it to be broken by NNN hoppings or 
otherwise. Ultimately, it is this broken chiral symmetry 
that is responsible for the sensitivity to selective disorder. 

We understand the topological transition as a transi¬ 
tion from the honeycomb (topologically non-trivial for 
finite Haldane mass) to the two uncoupled triangular 
(trivial) sublattices within the same symmetry class. The 
“localization transition” is more complicated but in the 
limiting cases it is also between a Chern insulator (topo- 
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logically non trivial) to a metal, again without changing 
symmetry or dimensionality. The localization transition 
in graphene with vacancies equally distributed among the 
two sublattices has been explored in some detail as a 
special example of chiral class 18 20 presenting some dif¬ 
ferences with the Anderson disorder. But the sensitivity 
of the localization transition to the sublattice selection 
has only been addressed more recently! 2 Reference [37 
discusses the topological transition in the case of the time 
reversal symplectic class (Kane-Mele model). Our find¬ 
ings are qualitatively similar in the equally diluted case 
but differ in the selected dilution case. Ref. [21] analy¬ 
ses the localization properties of a 2D chiral systems (in 
particular in graphene) without topological features with 
emphasis on the selected dilution case. Our results can 
not be directly compared with theirs because the NNN 
complex hoppings make a big difference. Apart from be¬ 
ing responsible for the topological non trivial features, 
in the case of selective vacancy disorder, NNN provide 
the connectivity of the unaffected triangular sublattice. 
In the extreme case when one sublattice is totally de¬ 
pleted, the remaining sublattice in our case is metallic 
with bands shown in Fig. [5] of Appendix |A| It will be in¬ 
teresting to analyze the scaling behavior of the Haldane 
model along the lines of Ref. [21] A recent calculation of 
the conductivity tensor in the disordered Haldane model 
suggests that an asymmetry between sublattices A and 
B can help to stabilize the Chern insulator,^ an intrigu¬ 
ing result compatible with our findings. It is interesting 
to note that the Haldane model with purely imaginary 
hoppings (cf) = 7r/2) has particle-hole symmetry so that 
it belongs to class D whose Chern number also belongs 
to We do not expect significant differences in the 
results for the case of a general value of <fi. 
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Appendix A: From honeycomb to two triangular 
lattices 

Dilution of a single sublattice leads, after 100% dilu¬ 
tion, to the triangular lattice. Since the latter has a single 
band, the Chern number must be zero. In this section 


we analyze whether the transition found for the selected 
dilution case can be attributed to the fact that, for some 
critical density of vacancies, the system is already more 
triangular lattice like and therefore the Chern number is 
zero. 

To address this possibility we study how the gap and 
Chern number evolves with the NN hopping t. For t = 0 
we end up with the triangular lattice. The spectrum and 
the eigenvectors are obtained by the 2x2 Hamiltonian 
given by 


where 




m( k) /(k) \ 
/*(k) ~m( k) ) 


m(k) = —2£ 2 [sin(k • ai) — sin(k • (ai — a 2 )) — sin(k • a 2 )], 
/(k) = £(l + e* k - ai +e ik ’ a2 ), 

and we have chosen ai = a(— 1,a/3)/2 and a 2 = 
a(l, \/3)/2 with a the lattice constant for the underlaying 
hexagonal lattice. The spectrum is easily obtained and 
reads 


£k = ±a/to( k) 2 + |/(k) | 2 . (AI) 

In the limit t/t 2 —>> 0 we get two decoupled triangu¬ 
lar lattices with imaginary hoppings related by com¬ 
plex conjugation, the spectrum being F k = d=ra(k). In 
Fig. [4] we show the spectrum F k = m(k) for the triangu¬ 
lar lattice with imaginary hoppings (left), the spectrum 
F k = ±m(k) for two decoupled triangular lattices with 
imaginary hoppings (center), and the spectrum given by 
Eq. when the system is perturbed by a finite t lead¬ 
ing to a finite gap (right). 

It is apparent that any finite NN hopping parameter t 
lifts the degeneracy that can be seen in the middle panels 
of Fig. [4j and opens a finite gap in the spectrum. Since 
we know by inspection that the gap remains open for any 
t/t 2 , including t/t 2 1 where the system in topologically 
non-trivial, then we conclude that it must be non-trivial 
for any finite t. This is explicitly shown in the left panel 
of Fig. [5j where we represent the Chern number for the 
Haldane model as a function of t/t 2 . The Chern num¬ 
ber has been obtained using Fukui’s method.^ So, the 
system is trivial only for t = 0, when the two triangular 
lattices are decoupled. 

Let us emphasize that we should not confuse the trian¬ 
gular lattice with homogeneous complex hoppings, which 
breaks time reversal symmetry, with the familiar triangu¬ 
lar (or hexagonal) lattice with real hoppings. The spec¬ 
trum for the latter reads 


F k = —t 2 (cos(k • a 2 ) + cos(k • ai) + cos[k • (ai - a 2 )]) 


= -h : 



+ cos (ak x ) 


Such energy spectrum is shown in Fig. [ 5 ] (middle and 
right panels). 
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Figure 4. Energy spectrum for the triangular lattice with imaginary hoppings (left), two decoupled triangular lat tices with 
imaginary hoppings related by complex conjugation (center), and coupled triangular lattices (right) as given by Eq. (A1). Top 
panels are the same as bottom panels along a particular path in the BZ. 



Figure 5. (left) Chern number for the Haldane model as a function of the parameter t/t^. (middle) Energy spectrum for the 
triangular lattice with real hoppings, (right) The same as the middle panel along a particular path in the BZ. 
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